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Motivated by recent experimental results, we study the effect of size reduction on half-doped man- 
ganite, Lao.sCao.sMnOa, using the combination of density functional theory (DFT) and dynamical 
mean field theory (DMFT). We find that upon size reduction, the charge-ordered antiferromagnetic 
phase, observed in bulk, to be destabilized, giving rise to the stability of a ferromagnetic metal- 
lic state. Our theoretical results, carried out on defect-free nanocluster in isolation, establish the 
structural changes that follow upon size reduction to be responsible for this. Our study further 
points out the effect of size reduction to be distinctively different from application of hydrostatic 
pressure. Interestingly, our DFT+DMFT study, additionally, reports the correlation-driven stability 
of charge-orbitally ordered state in bulk Lao.sCao.sMnOa, even in absence of long range magnetic 
order. 



PACS numbers: 68.65.-k,71.27.+a,71.15.Mb 

Size controls the physical properties of materials and 
can hence be employed to make materials functional. 
For strongly correlated materials, theoretical modelling 
of such phenomena is rare. Here, we pursue such a 
study taking the case of half-doped manganites. The 
charge and orbitally ordered state observed in half- 
doped manganites is one among the rich variety of fas- 
cinating phenomena exhibited by perovskite manganites 
Ri_ 2 ,A a; Mn03 (R= rare-earth element; A= alkali-earth 
element). ^ The charge-ordered (CO) state is associated 
with a real space ordering of Mn 3+ / Mn 4+ species in a 
1:1 pattern. It is accompanied by orbital ordering (00) 
and a structural change from orthorhombic to mono- 
clinic symmetry which gives rise to an insulating ground 
statePtSI The insulating CO state has been reported to be 
destabilized to a ferromagnetic (FM) metallic phase by 
various means that include a magnetic fieldP doping, bi- 
axial strain, pressure,^ and electric fieldP Very recently, 
it has been shown in a few experimental studies^HH] th a t 
the destabilization of the CO state can be achieved even 
through size reduction. This interesting phenomenon 
adds another dimension, namely size, to the problem. 
Size control is attractive from a technology point of view, 
which is achievable chemically in a low-cost way. The 
route through size control also opens up the possibility 
of exploring the tunability of the CO-00 state and the 
associated metal-insulator transition. 

In this letter, we study the effect of size reduction on 
the CO-00 state of Lao. 5 Cao.5Mn0 3 (LCMO) by us- 
ing a combination of First-principles DFT and DMFT. 
For the DFT calculations, we used a combination of two 
methods: (a) plane wave-based pseudopotentialspSl and 
(b) muffin-tin orbital (MTO) based on linear muffin-tin 
orbitaP (LMTO) and N-th order MTO (NMTO)PFor 
(a) we used projected augmented wave pseudopotentials 
with an energy cutoff of 450 eV. We used a spin polar- 
ized generalized gradient approximation (GGA)PSlFrom 
a self-consistent DFT calculation, a low-energy Mn-e s 
only model Hamiltonian was constructed using NMTO- 



downfolding techniqueP The corresponding Hubbard 
Hamiltonian defined in the downfolded NMTO basis was 
solved using DMFT, in the same spirit as previously car- 
ried out in in Ref. [T7| in the context of pure LaMn03. 
The low-energy model Hamiltonian consists of two e g or- 
bitals per Mn ion with intra-orbital Coulomb interaction 
U = 5eV and Hund's exchange J = 0.75 eV which are 
coupled by J = 1.35 eV to a (classical) spin representing 
the half-filled and immobile t 2g electrons. The DMFT 
equations were solved by Hirsch-Fye quantum Monte 
Carlo simulations^, and, because of the CO ordering, 
it was necessary to explicitly consider a site-dependent 
double counting correction P2a 

Bulk LCMO shows a CO transition at T co = 155 K. 
Below 155 K, the crystal structure is of monoclinic sym- 
metry (P2i/m) and an antiferromagnetic (AFM) order 
sets inP3 The magnetic structure, the so called "CE" 
order, consists of zig-zag FM chains that are coupled 
antiferromagnetically in the crystallographic ac plane. 
The ac planes are stacked antiferromagnetically along 
the crystallographic 6-direction. A noteworthy feature of 
the crystal structure is that while the Jahn- Teller (JT) 
distortion is sizable for the bridge-site Mn atoms (Mnl) 
with two long bonds along the FM chain and four short 
bonds, the corner-site Mn atoms (Mn2) on the zigzag 
chains have negligible distortion with nearly similar Mn- 
O bondlengths. Average Mn2-0 distance is smaller than 
that of Mnl-OPOur DFT calculations carried out on the 
experimentally measured structure^ henceforth, referred 
as S ea: , showed the CE insulating phase to be stable by 18 
meV/f.u. over the FM metallic solution. The calculated 
electronic structure in terms of density of states and mag- 
netic moments are found to be in good agreement with 
those reported previously in literature.^ 

In order to study the problem of nanoscale LCMO we 
first created a large supercell in the monoclinic structure, 
from which a cluster of diameter 2 -3 nm having approx- 
imate spherical shape was cut out (cf Fig. 1, left panel), 
in which we carried out a full structural optimization. 
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a = 5.47, b = 7.58 
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FIG. 1: (Color online) Left Panel: Constructed nanocluster of 
LCMO. In magenta (dark grey), brown (light grey) and blue 
(deep, dark grey) we show the MnOg octahedra belonging to 
the outer most surface layer, next to the surface and the core 
region. Right Panel: The structural unit, highlighted in green 
(grey), chosen out of S na no, used for building up of Smodei- 
The big and small spheres in the unshaded region represent 
Mn and O atoms respectively. 



Following this procedure, the 2 nm cluster contains a to- 
tal of 370 atoms and the 3 nm cluster contains a total of 
700 atoms, pushing it to a limit of our DFT structural 
optimization. In the construction of the clusters, care 
has been taken to maintain the stoichiometry as closely 
as possible. For the cluster calculation, a simple cubic 
supercell was used with periodic boundary conditions, 
where two neighboring clusters were kept separated by 10 
A, which essentially makes the interaction between clus- 
ter images negligible. The positions of the atoms were 
relaxed towards equilibrium, using the conjugate gradi- 
ent technique until the Hellmann-Feynman forces became 
less than 0.001 eV/A. 

The considered DFT cluster sizes are smaller than the 
experimental realizations^ of sizes 15 nm. Hence only 
the inner region of the above constructed clusters of 2-3 
nm size, is expected to mimic the prototypical behavior 
of the experimentally studied clusters. In order to under- 
stand the consequences of the size-controlled structural 
changes for such relatively larger clusters, we hence con- 
structed a model bulk system, which we refer as S mo dei- 
It is built out of the structural units belonging to the 
innermost core and the next to the core layer of the op- 
timized LCMO in the nanoscale geometry, (referred as 
Snano) 7 as shown in the right panel of Fig. 1, and sub- 
sequently imposing the symmetry considerations. The 
detailed procedure of construction of model system is 
explained in the supplementary information (SI). Con- 
struction of Smodei loads to consideration of the the local 
oxygen environments around Mn atoms as well as the 
tilt and rotation connecting two MnOg octahedra, same 
as that in core region of S„ a „ . The lattice parameters 
and the Mn-0 bond lengths of S m odei are compared to 
the bulk structure in Table I. The detailed structural in- 
formation can be obtained in SI. Note, for the bulk, we 
have considered the theoretically optimized structure, re- 



TABLE I: Lattice parameters, and Mn-0 bondlengths (in 
A) of Smodei in comparison to Sbuik- The entries for Mn- 
O bondlength from left to right correspond to that along the 
FM chain, between the FM chains, along 6-direction and the 
ave rag e. Mnl atoms are of nominal valence 3+ and Mn2 of 

4+ pr 



ferred as Sbuik, hi order to compare with parameters of 
Smodei m the same footing. We find that the lattice pa- 
rameters of S mo dei show substantial reduction compared 
to those of the bulk system. The change in a param- 
eter appears to be the largest with a change of about 
0.20 A, with moderate changes in b and c parameters, 
of 0.09 A's. Qualitatively, this trend of reduction in lat- 
tice parameters and also the nature of reduction agrees 
very well with the crystal structure data extracted from 
X-ray diffraction of nanoclusters of LCMO of 15 nm size 
[cf. Fig 4 in ReffTUTa)]. We note that the reduction in 
lattice parameters in the model structure gave rise to 
about 6% reduction in the volume compared to that of 
the bulk system; the first experiments^ report a 2% re- 
duction. The 6% reduction was obtained for S mo d e i con- 
structed out of S nano of 3 nm, while a similar procedure 
for S nano of 2 nm, gave rise to larger volume reduction 
(~ 8 %). This indicates the volume reduction to increase 
upon decreasing cluster size, justifying the difference be- 
tween the obtained volume reduction on 2-3 nm cluster 
and experimentally observed volume reduction on 15 nm 
cluster. One of the important structural quantities is 
the orthorhombic strain: OS n = 2 \ c , a \ gives the strain 

II (c+a) o 

in the ac plane, while OS i = 2 ( a+c ~v^>) j s that along 

(a+c+V2b) 

the 6-axis. For Sbuik, the orthorhombic strain is highly 
anisotropic with a negligible value of OSn (~ 0.002) and a 
high value of OSj_ (~ 0.021). For S mo dei, we nn d instead 
the orthorhombic strains to be comparable (OSn ~ 0.02 
and OSj^ s» 0.01). This trend is also in very good agree- 
ment with experimental results.^ It confirms that our 
constructed model structure captures the essential struc- 
tural changes in the nanoscale surprisingly well. This 
proves that the role of surface beyond what is already 
taken into account in construction of model structure, is 
small. 

Next, we calculated the electronic structure for S m odei 
and compared it with that of Sbuik- For understand- 
ing the results, let us note that the difference between 
the average Mnl-0 and Mn2-0 bondlengths is smaller 
in S mo dei than in Sbuik- This leads to the expectation 
that the charge disproportionation (CD) between Mnl 
and Mn2 sites to decrease in S mo dei- Furthermore, we 
note that for Sbuik, the difference between the longest 
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FIG. 2: (Color online) DFT DOS, projected onto Mnl-d 
(black solid line) and Mn2-d (green/light grey line) states 
calculated for the CE phase of Sbuik (left panel) and the FM 
phase of S mo d e i (right panel). The zero of energy is set at Ef- 



and the shortest Mn-0 bondlengths is large for Mnl and 
tiny for Mn2. This gives rise to the crystal field splitting 
(A) between the two Mn-e 9 states, Mn-3z 2 — r 2 and Mn- 
x 2 — y 2 , as large as 0.63 eV for Mnl sites and less than 
0.02 eV for the Mn2 sites. In contrast for S mo( z e z, the 
bond length differences are much more similar for both 
types of Mn sites. This is reflected in similar A's for the 
nano model, i.e., 0.15 eV for Mnl sites and 0.10 eV for 
Mn2 sites. Together these two effects weaken CO as well 
as OO in S m odei- This ordering is important to stabi- 
lize the AFM structure found for the bulk. With charge 
and orbital ordering weakened, we find instead FM to be 
stable by 20 meV in S mo dei ■ This r esult is in accordance 
with the experimental observations!^^ The microscopic 
origin of the size controlled transition from AFM to FM, 
therefore, can be traced back to the size-induced struc- 
tural changes. 

Fig. [2] shows the density of states (DOS) of Sbuik with 
AFM ordering of Mn spins, in comparison to that of 
Smodei with FM ordering. Considering the DOS for Sbuik, 
the crystal field splitting between Mn-3z 2 — r 2 and Mn- 
x 2 — y 2 is clearly seen. In the majority spin channel, Mn- 
3z 2 — r 2 states at Mnl site are more occupied than the 
Mn-3z 2 — r 2 states at Mn2 site, giving rise to CD between 
Mnl and Mn2. We also find OO at Mnl sites with a pref- 
erential occupation of Mn-3z 2 — r 2 over Mn-x 2 — y 2 . The 
CO, although incomplete, together with the AFM spin 
ordering gives rise to an insulating solution with a small 
but finite gap at F,p already at the DFT level. Consider- 
ing the DOS of Smodei, we find that the splitting between 
Mn-3z 2 — r 2 and Mn-1 2 — y 2 is less pronounced and the 
Mnl-(i and Mn2-d states to be similar. The reduced A 
together with the increased bandwidth, compared to the 
bulk structure, drives Smodei to metallicity with a finite 
density of states at Fermi energy (Ep). The increased 
bandwidth is caused by the reduction in volume as well 
as by the FM ordering which allows hopping processes 
within a double exchange model!^ 

In order to take into account the influence of the miss- 
ing electronic correlation in GGA, we did paramagnetic 
DMFT calculations for Sbuik and S m odei structures, con- 





Bulk 


Nano-model 


Mnl(l) 
Mnl (2) 
Mn2 


0.87 (0.50) 0.01 (0.11) 
0.85 (0.47) 0.01 (0.12) 
0.04 (0.15) 0.09 (0.25) 


0.52 (0.31) 0.09 (0.20) 
0.72 (0.38) 0.04 (0.19) 
0.16(0.21) 0.16 (0.25) 



TABLE II: Orbital occupancies for Mn-3z 2 - r 2 (first en- 
try) and x 2 — y 2 (second entry) states calculated within 
DFT+DMFT for different inequivalent classes of Mn atoms in 
the unit cell for Sbuik and S mo dei- In brackets we give the cor- 
responding occupancies for the one particle, low-energy DFT 
Hamiltonian. 
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FIG. 3: (Color online) DFT+DMFT spectral densities for e 9 
states of three inequivalent classes of Mn (solid lines), cor- 
responding to Sbuik and S mo£ j e ;. The dashed lines represent 
the corresponding DFT DOS. The lines colored as black and 
green (light grey) correspond to 3z 2 — r 2 and x 2 — y 2 states, 
respectively. 



sidering the low-energy Mn-e 9 only Hamiltonian derived 
out of DFT. Table II lists the the orbital occupations of 
the three types of inequivalent Mn sites. For Sbuik, a l~ 
ready in the DFT (in brackets) the two inequivalent Mnl 
( "Mn 3+ -like" ) sites (Mnl(l) and Mnl(2))_are more oc- 
cupied than the Mn2 ("Mn 4 +-like") sitesPl Besides the 
CD, there is, as mentioned above, a DFT orbital order. 
Electronic correlations enhance both kinds of ordering 
dramatically, making CO and OO nearly complete. This, 
establishes the correlation driven stability of CO and OO 
with almost complete CD in a paramagnetic phase. The 
almost complete CD and enhanced orbital polarization 
(OP) at the Mnl sites results in a gap at the chemical 
potential in the DFT+DMFT spectral function for the 
bulk structure, as shown in the left panel of Fig. [3j even 
without spin ordering. Note the opening of charge gap 
is stabilized by long-range Coulomb interactions which, 
within DMFT, reduce to their Hartree contribution. This 
effect is taken into account in our DFT+DMFT calcula- 
tion on the GGA level. Compared to the DFT spectra, 
spectral weight is transferred to high frequencies in the 
form of Hubbard bands, opening a gap at the chemical 
potential. For the DFT calculation on the other hand, 
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the charge disproportionation is incomplete, and the in- 
sulating solution is obtained only by assuming the AFM 
spin ordering. Turning to S mo deU the DFT occupan- 
cies show Mn 3+ -like and Mn 4+ -like sites to be similar 
with only a weak CD. The inclusion of correlation effect 
through DMFT enhances CD to some extent following 
the trend seen for Sbuik- However, CD remains incom- 
plete with an average occupation of Mn 3+ -like and Mn 4+ - 
like sites of 0.6-0.7 and 0.3 respectively, in comparison to 
0.9 and 0.1 respectively, obtained for Sbulk- This conclu- 
sively establishes that size reduction leads to weakening of 
charge disproportionation. This in turn leads to metal- 
lic DFT+DMFT solution for S mode i with finite weight at 
the chemical potential, as shown in right panel of Fig. 3. 
Note, although S nono does not maintain the strict stoi- 
chiometry, the constructed S mo dei is strictly stoichiomet- 
ric, pointing the fact that destabilization of CO is driven 
by the structural changes due to size confinement, rather 
than due to deviation from half-doping. 

As one of the major structural changes upon size re- 
duction is the volume compression, it is worthwhile to 
compare the structural and electronic changes induced by 
size reduction to those occurring under hydrostatic pres- 
sure. To this end, we carried out calculations of LCMO, 
with uniformly reduced lattice parameters with 6% re- 
duced volume, the atomic positions being optimized in 
DFT, referred as structure S press . The details of the 
optimized structure is given in SI. Following the self- 
consistent DFT calculations on S press , the Mn-e g only 
low-energy Hamiltonian was constructed and the corre- 
sponding Hubbard Hamiltonian was solved using DMFT. 
Compared to S mo dei , first of all, we find that at the DFT 
level, CD and OP is much weaker, even though the vol- 
ume is the same. With this less polarized starting point, 
all Mn sites are filled with « 0.5 electrons. In this situa- 
tion, electronic correlations are less relevant. The DMFT 
orbital and site occupations remain very similar to the 
DFT values with 0.4 - 0.6 electrons/site, and the system 
is far away from a metal-insulator transition (MIT) . This 
leads us to conclude that the nanoscopic system is much 
closer to a MIT than bulk Lao.sCao.sMnOs under hydro- 
static pressure.^ The size reduction and application of 



hydrostatic pressure, therefore, should be considered as 
two very different routes. 

In conclusion, using DFT calculations combined with 
DMFT, we have studied the effect of size reduction on 
charge-orbital order in half-doped LCMO manganites. 
Our study indicates that the size reduction leads to 
substantial reduction in volume as well as a change in 
the nature of the orthorhombic strain. The structural 
changes under size reduction lead to a weakening of both 
charge and orbital ordering, making the ferromagnetic 
metallic state energetically favorable compared to the 
"CE" type antiferromagnetic insulating state, which is 
the ground state of the bulk structure. While such effect 
has been observed, the experimental situation is faced 
with difficulties, like possible presence of impure phases, 
the grain boundaries, non-stoichiomctry. Our theoret- 
ical calculations were carried out considering nanoclus- 
ter in isolation, and therefore, devoid of such complica- 
tions. Through construction of model structure, the is- 
sue of non-stoichiometry was avoided. Furthermore, the 
effect of size reduction turned out to be very different 
from that of pure hydrostatic pressure. We predict the 
nanoscopic system to be close to the MIT in compari- 
son to the system under hydrostatic pressure with the 
same amount of volume reduction. Increasing the size of 
the nanocluster, one would expect to drive the system 
closer and closer to MIT. It would therefore be possi- 
ble to tune the LCMO system to the verge of MIT, and 
thereby, achieve a large magnetoresistive response un- 
der small magnetic fields. Finally, while we carried out 
our investigation on LCMO, the destabilization has been 
predicted for Pro.sCao.sMnOa toop^ hinting to oberse- 
ved effect to be a more general one. This will be taken 
up in a later study. 
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